UCLA/01/TEP/15 

Vortex waistlines and long range fluctuations 



Tamas G. KovacsQQ 

NIC/DESY Zeuthen, Platanenallee 6, D-15738, Zeuthen, Germany 

e-mail: kovacs@ifh.de 

and 

E. T. TomboulisFj 

Department of Physics, UCLA, Los Angeles, CA 90095-1547 
e-mail: tombouli@physics.ucla.edu 



Abstract 



We examine the manner in which a linear potential results from fluctuations due 
to vortices linked with the Wilson loop. Our discussion is based on exact relations and 
inequalities between the Wilson loop and the vortex and electric flux order parameters. 
We show that, contrary to the customary naive picture, only vortex fluctuations of 
thickness of the order of the spatial linear size of the loop are capable of producing a 
strictly linear potential. An effective theory of these long range fluctuations emerges 
naturally in the form of a strongly coupled Z(N) lattice gauge theory. We also point out 
that dynamical fermions introduced in this medium undergo chiral symmetry breaking. 
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1 Introduction 



Center vortices are widely believed to be the most important degrees of freedom for 
confinement in Yang Mills theories. (For a recent review see Attempts to isolate 
and compute the vortex content in the path integral at large j3 have been a very active 
area of study in the last two years. (For a general view of the various approaches and 
issues involved, see 0). 

A popular plausibility argument for confinement by vortices goes as follows. To 
estimate the vortex contribution to the expectation of a given Wilson loop, one assumes 
that typically many thin vortices — i.e. thin compared to the physically large (3> lfm 2 ) 
Wilson loop area A — link with the Wilson loop. A simple estimate is obtained by 
subdividing the surface spanned by the loop into n = A/a pieces of equal area a and 
assuming that a vortex piercing a given small area piece is present with probability p. 
Provided that the vortices linking with the loop at different places are independent, 
the Wilson loop expectation is then given by 



which is the desired area-law. 

In the present paper we point out that while this simple picture indeed gives the 
correct area law, it is fundamentally incorrect because it is not supported by the un- 
derlying full non-Abelian gauge theory. The response of the vortex distribution in the 
vacuum to the introduction of the external probe represented by a Wilson loop should 
be such as to minimize the effective quantum action (free energy) of the system. To 
examine this we make the above arguments more precise by introducing a set of vortex 
containers linking with the Wilson loop in the spirit of ||. Any set of vortex contain- 
ers can then be used to obtain an upper bound to the Wilson loop by making use of 
rigorous inequalities mostly based on reflection positivity. We show that thick vortex 
containers of thickness of the order of the (shorter) side length of the Wilson loop yield 
the 'best' and in fact the only strict area- law upper bound. This shows that vortices do 
not have some fixed characteristic thickness but the most important vortex fluctuations 
disordering a Wilson loop of a given size have a core thickness of the order of the linear 
size of the Wilson loop. It also means that generally vortices cannot be unambiguously 
located as individual objects in every given single gauge field configuration. The only 
meaningful quantity is the number of vortices (mod N) linking with a given Wilson 
loop, and this is a well-defined gauge invariant quantity. In the remainder of the paper 
we show how the physical picture implied by the previous discussion leads to a simple 
effective Z(N) gauge theory for the long distance center degrees of freedom. We then 
observe that the introduction of dynamical fermions in this effective strongly coupled 
system results in chiral symmetry breaking. 

The paper is organized as follows. In Section 2 we review and extend various 
relations and inequalities relating the electric free energy (Fourier transform of the 
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vortex free energy) to the Wilson loop. These relations form the basis for our discussion 
of vortices and the emergence of a linear potential in Section 3. Section 4 further 
discusses the physical picture and introduces the effective theory. Section 5 contains 
some conclusions. 

2 Electric flux inequalities 

After a brief description of our notation, in this Section we introduce the basic electric 
flux inequalities. These are already interesting in their own right since they yield a 
rigorous upper bound on the Wilson loop in terms of the electric flux free energy. 

We work on a hypercubic lattice A of length L M in spacetime direction fi = 1, . . . ,d. 
We assume the standard Wilson formulation of lattice gauge theory with SU(2) group- 
valued link variables, the Wilson action and periodic boundary conditions in all direc- 
tions. Expectation values of observables are defined as 



where the integration is over all the group- valued link variables and Z, the partition 
function, is the same integral without the operator insertion O. 

Let us denote by C[V Miy ] the operator that flips the sign of the coupling (introduces a 
Z{2) 'twist') on a coclosed (d— 2)-dimensional set of plaquettes V^ v winding around the 
periodic lattice in the directions perpendicular to the /^-directions. The expectation 
value of this operator defines the vortex free energy: 



The twist amounts to a discontinuous gauge transformation with multivaluedness in 
Z{2), i.e. forces the presence of a tti(SU(2)/Z(2)) vortex wrapped around the periodic 
lattice. As indicated by the notation, the expectation depends only on the directions in 
which V winds through the lattice, not the exact shape or location of V. This expresses 
the mod 2 conservation of flux. Indeed, the twist —1 on the plaquettes forming V can 
be moved to the plaquettes forming any other homologous coclosed set V ' by the 
change of variables U\, — > — Uf, in the numerator in (^) for each bond b in a set of bonds 
cobounded by V U V '. By the same token (|3]) is invariant under changes mod 2 in the 
number of homologous twisted coclosed sets introduced in A. A simple consequence of 
this is that 



We will assume that, for sufficiently large and dimension d < 4, the vortex 

free energy (|3|) behaves as 




(2) 




(3) 



(o[v^)o[v; v ]) = i. 



(4) 



i^H~(Il L x ) exp(-p(/3)|Vl) 



(5) 
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where {A^l = L^Ly. This is the optimal behavior under exponential transverse spread- 
ing (creation of mass gap) of the flux introduced by the twist on V, with p approaching, 
at least asymptotically, the exact linear potential string tension. This behavior is ex- 
pected by physical reasoning and explicitly seen in the strong coupling expansion 
||. Recently, it became possible to demonstrate this in numerical simulations at large 
(3 ||, f7|. The behavior (^) for a 'vortex in a box' is essential for our argument in the 
following. 

The Z(2) Fourier transform of (Q), 

exp(-F d ) = (i(l-0[V^])), (6) 

gives the corresponding dual (w.r.t. the center) color-electric free energy. The mod 2 
conservation of the magnetic flux is now expressed by the projection property 

{\{i-o[v» v \)\{i-o[v; v \)) = -o[v^)), (?) 

as is easily seen by using Eq. (|). 

Consider now a rectangular Wilson loop C placed, say, in the [12]-plane. Let V be 
a coclosed stack of plaquettes winding around the periodic lattice in the perpendicular 
\i — 3, . . . , d directions and through C, and insert unity in the numerator in the Wilson 
loop expectation W[C] = ( |tr[/[C]y in the form 

1 = \{1-0[V])+ 1 -{1 + 0[V]). 

Then 

W[C\ = {\tTU[C}\(l-0[V])) + (ltTU[C}\(l + 0[V})) 

= ( \t.U\C} ~( 1 - 0[V] ) ) + ( \%rU[C] i( 1 - 0[V] ) ) (8) 

Here V ' is another coclosed stack of plaquettes winding around the lattice in the per- 
pendicular [A — 3, . . . , d directions but not threading through C . The second equality 
in @ is obtained by making the change of variables Ub — > —Ub in the second term in 
the first equality in (||) for each bond b in a set of bonds cobounbed by V and V '. This 
'moves' the twist (-1) on the plaquettes forming V to those forming V '. This set of 
bonds necessarily involves one (or an odd number of) bond(s) on C, which results in 
the minus sign in the second term in the second equality. (^) is represented graphically 
as 




(9) 



where a filled square stands for the operator |( 1 — 0\V] ), with V crossing the two- 
dimensional plane containing the loop at the location of the filled square, and winding 
around the lattice in the remaining d — 2 perpendicular directions. 
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Simple identities like (^), or ([LED , QlTD below, serve as the starting point for deriving 
relations between the Wilson loop and vortex free energies by use of reflection positiv- 
ity. Given a reflection r about a (d — l)-dimensional lattice hyperplane, one defines 
an antilinear mapping 9 on functions of the bond variables by #.F[{{7&}] = _F*[{[/ r (,}]. 
By the reflection positivity properties of the LGT action |J, this induces a positive 
semidefinite inner product on the space of configurations allowing the use of the cor- 
responding Schwarz inequalities. Thus, starting with (P), consider a reflection about a 
d— 1- dimensional hyperplane 7r perpendicular to the 'vertical' loop legs and containing, 
say, the top 'horizontal' leg of the loop C. One then has the inequality 




(10) 

for the first term, and a similar statement for the second term in @. In fllTf ) we made 
use of (0)0 and |trl = 1. The loop has now doubled in size along the direction of 
one of its legs. Proceeding now by repeated reflections in hyperplanes containing one 
of the legs of the loop resulting from the previous reflection, and use of (|7|), one may 
eventually completely eliminate all tr^ factors from the loop operator by virtue of the 
lattice periodicity. Applying this procedure to both terms on the rhs. of Eq. ([|) one 
obtains 

/ \ \A C \/\A± 2 \ 



W[C\) <2^exp(-F cl )j , (11) 

where \Ac\ is the minimal area bounded by the Wilson loop. If then the vortex free 
energy behaves as in ([5]), (|Il]) implies area law for the Wilson loop. 

Note that the result (|TTD manifestly incorporates mod 2 conservation since of course 
(P) and (HD do. This is an important point that we now explore a bit further. Any 
multiple factors of |( 1 — 0[V] ) occurring in the derivation above were eliminated by 
(0). Suppose instead that we keep a number of such factors to make contact with the 
naive picture of a Wilson loop pierced by several independent vortices. So imagine that 
we subdivide the 2-dimensional plane A containing C into large squares of side length 
/, i.e. we view A from a coarse lattice of spacing /. We denote A by A' when viewed 
from the coarse lattice. Rewrite the identity (§) using (0) in the equivalent form 

W[C] = ( ~tr[/[q n \{ 1 - )i > + ( \toU[C\ II ~( 1 - 0[V/} ), ) (12) 

i 3 

where the product in the first term includes one factor (indexed by i) for every large 
square in A ( every plaquette in A') tiling C (figure [ft), and the identical arrangement 

4 Notice that (Q) remains valid in the presence of additional operators in the expectation, as long 
as moving V ' to V by a change of variables does not affect the additional operator. In particular, this 
is the case for a Wilson loop when neither or both V and V ' link with the loop. 
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Figure 1: Arrangement for reflections in first term in (|T|) to obtain (|l"3"D (see text). 



translated outside C in the second term. One now applies reflection positivity to 
repeatedly reflect about (d — l)-dimensional hyperplanes ir[, or tc' 2 , perpendicular to 
the /i = 1, or 2-direction, respectively, containing bonds in A' (figure |l|). In this way 
one derives 

W[C]<2U ( I] ^(1-0[V,]),) 1/|A '' (13) 

* keA' 

where in the product inside the expectation there is one factor for each plaquette in 



A'. (fyf) now has the form of a typical 'chessboard estimate' inequality (see e.g. |§, 
0). It is of course equivalent to (|TT| ) by (|7|) since \A'\ = \A\/l 2 , and the number of 
factors in the outside product equals \A' C \ = \Ac\/l 2 - 

Returning to (§) and noting that 



\%rU[C]\{l-0[V] 



UxU[C]\{l-0[V '} 



\trU[C] U 1 - 0[V] ) M 1 + 0[V '])) + ( \trU[C} |( 1 + 0[V] ) |( 1 - 0[V ' 



\tvU[C]\(l-0[V]0[V '] 
\trU[C]\{l-0[V"])) 
one obtains the alternative identities represented graphically by 



(14) 
(15) 



V 



V 



W[C] 



,x d 



(16) 
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(|16|) is the graphical representation of (fH]). (p~5|) , graphically depicted in fllTD , is then 
obtained from flTil) by simply merging V and V ' in the coclosed set V " linking with 
C by another shift in integration variables and mod 2 flux conservation.^ 

Alternatively, (|i~6l) , (|17]) may be directly obtained as follows. Insert 1 = J*v 2 ) ^7; 
where 7 G Z(2), in the numerator in the expectation and make a shift of inte- 

gration variables Uf, — > 7 [/& for all b G B, where B is a set of bonds whose coboundary 
is V U V ' or V ". The result is fllED , or (0), respectively. 



3 The Wilson loop in terms of vortex containers 

In this Section we introduce a set of vortex containers linking with a given Wilson 
loop. Using the relations derived in the previous Section, we obtain an upper bound 
on the Wilson loop in terms of vortex fluctuations occurring in the given set of vortex 
containers. We then ask the question, which set of vortex containers gives the "best" 
(i.e. the lowest) upper bound to the Wilson loop. It turns out that the favored set of 
vortex containers has only one single thick container utilizing an area of the order of 
the area spanned by the Wilson loop. This yields a strictly linear lower bound to the 
heavy quark potential, whereas a collection of many thin vortex containers results in 
a suppression of the potential by a logarithmic factor. 

By (|7|) and footnote f|, one can insert multiple §( 1 - 0[V] ) factors in QT3p , flT7p 
corresponding to a collection of coclosed sets {Vj} linking with the loop C. Imagine 
enclosing each V« in a 'vortex container' 0, i.e. a sublattice Aj C A containing Vj 
and wrapping around C (figure 0). Imagine integrating over the bond variables in the 
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Figure 2: Vortex container Aj enclosing coclosed set of plaquettes Vj (shaded) linking 
with Wilson loop C. d = 3, or 3-dimensional section in d = 4. 



5 V U V ' ~ V " is topologically trivial w.r.t. the lattice T d topology, but not w.r.t. the obstruction 
of the loop C. 
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interior of each container, keeping bond variables on its boundary <9Aj fixed, and define 

1( 1 zffflW " 



where 



z£\u aAi )=[ II ^ II exp(§( T l) v ^tr[/ p ) (19) 

with the characteristic function Vj[p] = 1 if p e V i; otherwise. ^^(^aJ is of course 
simply the partition function for (the interior of) Aj. Then one obtains |3|]: 

W[C\ = (\trU[C\ UhAUaAj)- (20) 

i 

( |13D can be rewritten in the same way: 

W[C}<2U ( J] f l A k (U 9Ak )) 1/lA \ (21) 
» fceA' 

where now each container containing the coclosed set 14 is of transverse area I 2 and 
wraps around the lattice in the longitudinal directions. 

Note that /a^^aJ, (|18|) , is nothing but (fj|) (Z(2) FT of vortex free energy) now 
defined on a lattice Aj (the vortex container) with fixed (instead of periodic) b.c. in the 
transverse directions, but still periodic b.c. in the longitudinal directions (in which the 
Aj's have torus topology by construction). Vortex containers including integration over 
fluctuations (summation of entropy effects) in their interior are a convenient device for 
discussing scales larger than their thickness in terms of free energy costs rather than 
the action of individual configurations. 

From ( pOl) one can now, trivially, obtain the bound: 

W[C]<l[m^\f Ai (U dAi )\, (22) 

where the maximum is taken over all values of the bond variables on the boundary 
<9A,. 



makes a direct connection with the naive notion of independent vortices wind- 
ing through a large Wilson loop, resulting in disorder, and area law, provided they are 
able to grow sufficiently long to pierce through the loop at any point. For this to be 
possible the vortices must be allowed to grow sufficiently thick to keep their free energy 
cost fixed as their length increases with that of the Wilson loop. Let di be the size 
of Aj in each of the two directions transverse to the set Vj used in its definition; its 
longitudinal size is given by |Vj|. We must first assume that all the vortex containers 
are thick enough to reach the regime where 

-H^(U dAi )/z A +\u dAi )) ~ |V 4 |exp(-p< 2 ) (23) 
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for all Ugji,- So to keep the free energy cost of each vortex less than a fixed value /, 
we need 

d? > ~ ^ (Pfj ■ (24) 

Then also max f^iUdhi) < kO- ~ e ~*) < cons t- But \Vi\ (in d dimensions) is of the 
order R d ~ 2 for linkage through points away from the perimeter of a rectangular loop 
of side lengths T and R, T > R. Such a loop can then accommodate ~ RT / In R 
containers wrapped around it. Thus (p2|) gives a confining but not quite purely linear 
potential 

V(R) > const RJ In R. (25) 

The same reasoning, and consequent failure to produce a purely linear potential, 
applies to the familiar argument for confinement by vortices, outlined in the introduc- 
tion, which tacitly underlies or is implied by many discussions in the literature. One 
assumes randomly distributed vortices of a certain thickness and basically arbitrary 
length. It is crucial that one assumes that the cost and hence the probability for vor- 
tices to link anywhere with the loop is fixed for any Wilson loop size (cp. ([[])). One 
then considers one vortex linked with a large Wilson loop. With the vortex thickness 
assumed much less than the loop's linear dimensions, one now sums over all positions 
of intersection with a surface spanning the loop. With the above assumptions, this 
produces a factor proportional to the loop area. One then sums over all intersection 
points for two independent vortices linking with the loop, and so on. This clearly 
exponentiates generating area law: 

W[C] ~ 1 + (-1) k\A c \ + (( ~ 1) ^ cl)2 + ■ ■ • = exp( -k\A c \ ) • (26) 

We now see that purely linear confinement is obtained this way only by adopting a 
non-interacting gas picture, and ignoring the actual free energy requirements for having 
vortices of sufficient length link anywhere with a large loop: the type of discussion just 
given above for fl22D applies to each term in such a summation. Thus, for one vortex 
linking with the large loop of side lengths T and R (T 3> R) , a vortex cross section area 
of order In R is required; otherwise, linking anywhere far away from the perimeter for 
fixed, bounded vortex free energy cost /, as required by the argument, is not possible. 
This leads at best to (|25|) , not (p(f) . 

The problem arises because one treats the vortices as localized and independent. 
For sufficiently thick vortices free energy costs are indeed correctly estimated in mag- 
nitude as above, i.e. fl23|). Thus, if one imagines each vortex enclosed in a vortex 
container of fixed, but sufficiently large, width d, the exponential transverse spreading 
~ exp(— pd 2 ) renders the overall vortex bulk free energy cost inside insensitive to the 
exact values of the gauge fields on the container boundary. The vortex, however, is 
surrounded by the pure gauge long tail that encodes its nontrivial topology, and flux 
quantization. This tail incurs no additional action cost, but is of infinite range and 
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communicates the presence of nontrivial topological flux inside the container to every- 
where outside. So even though the gauge field values on a thick container's boundaries 
are irrelevant for estimating the bulk free energy cost inside, they are very much rele- 
vant for signalling the presence of a vortex inside to other vortices or other topological 
obstructions outside. 

This acts like an 'irreducible' interaction between vortices that acts at all distances, 
and enforces flux conservation mod N. This interaction allows a system of vortex 
excitations to adjust the amount of flux spreading, i.e adjust the thickness of vortices 
to minimize the free energy of the system. The thickness of vortex cores then is not 
fixed, but is adjusted relative to their length as required by the presence of other 
vortices and/or other obstructions (e.g. Wilson loop legs) sensitive to the presence 
of topological Z(N) flux (figure |3|). This means that in general vortices cannot be 
considered isolated, and a definite number of vortices, specified more precisely than 
mod N, cannot necessarily be unambiguously assigned to every configuration. 




Figure 3: Vortices linking with Wilson loop C. Cores shown in darker shading, and long 
range (infinite extent) pure gauge tails in light shading: (a) Well-separated vortices of 
approximately fixed width; (b) Configurations of equivalent mod 2 flux having one 
thickened vortex lowering potential energy; (c) Configurations contributing essentially 
as in (b) showing the possible ambiguities in counting vortices beyond mod N. 

Thus, in the presence of the Wilson loop source, the optimal configurations for 
the system are not those of multiple isolated linked vortices, each of some fixed free 
energy cost (figure |(a)), hence length |V| ~ R (d ~ 2) and fixed width d 2 ~ lni? < R,T. 
It is more advantageous, in terms of free energy cost, for multiple linking vortices to 
thicken and merge, the total topological flux being conserved mod N (figure H(b)). 
Since the Wilson loop operator is affected by the topological flux through it only mod 
N, this should optimize the expectation. But then the picture and expansion in terms 
of groups of isolated independent vortices (|26| ) is no longer applicable. 
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Similarly, (p2|) leads to (|25| ) because it is obtained by assuming the vortex inside 
each container as completely isolated and independent of all the others. The exact 
expression (|20| ) holds for any number of factors in the product inside the expectation 
(reflecting mod N conservation). In view of the above discussion, one may as well 
combine containers into ones as thick as possible by integrating over the boundary 
fields of neighboring containers. With T ^> R, this amounts to taking containers in 
the product in (p0|), (p2|) having transverse area ~ R 2 , and longitudinal extension 
~ (const R)( d ~ 2 \ now gives 

const 

V(R) > const R — (\nR + const) (27) 

R 

replacing (|25|) . For loops with T ~ R, basically only one vortex container is needed, 
which means that strict linear potential arises essentially from thick vortex fluctua- 
tions spanning the entire loop area. It is interesting to note that such thick vortices 
also produce nonperturbative 1/R contributions (at scales outside the short distance 
perturbation theory regime). 

Inequality ( p2]) is actually rather crude. The inequalities (|ll|), (0) following from 
reflection positivity are much more powerful because of the exponents that allow esti- 
mates uniform in the lattice size. They give directly pure area law. As is easily seen, 
this is true even if one further crudely bounds the r.h.s. of ( pT3[) from above as done in 
(f2"2]), since InL^/A' — > as |A| — > oo. There is no real reason for doing this though.^ 



4 Long range vortex fluctuations — effective theory 

The above discussion, based on exact relations and inequalities between the Wilson 
loop and the free energy order parameters, indicates that an effective picture of the 
long distance confining fluctuations as isolated, independent vortices winding over long 
distances - in other words, as some kind of an approximately dilute or weakly interact- 
ing vortex gas - is not generally applicable. It does not take properly into account the 
relation between length and thickness of a vortex implied by the cost in free energy 
necessary to create the vortex in the first place, nor the correlations between vortices 
caused by their long range (topologically nontrivial) pure gauge tails. These correla- 
tions are present irrespective of the separation and enforce the mod N conservation of 
topological flux. Even though the cost diminishes exponentially with the transverse 
thickness of a vortex (creation of mass gap), these effects must still be properly ac- 
counted for if vortices of basically arbitrary length are to be present in the vacuum. 
These effects then generally tend to cause neighboring thick vortices to thicken further 
and merge since this lowers the free-energy cost for the vortices (above any background 

6 (|l^), however, might serve as a starting point for more sophisticated approximations where the 
behavior (|5|) for the vortex free energy under spreading of flux is at least partially derived rather than 
assumed. 
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of fluctuations that may be present) for the same mod N total flux in the system. Thus, 
our previous discussion implies that the long distance linear potential should not prop- 
erly be viewed as arising from the fluctuations caused by a gas of independent vortices 
winding through the loop, each of thickness much less than the linear dimensions of 
the loop. Rather, the fluctuation is more accurately described as that due to vortices 
of thickness comparable to the (shorter) linear loop dimension. Negative values of the 
loop occur then with almost exactly equal measure weight as positive values (no vortex 
(mod 2)). This optimizes the expectation to be as small as possible, i.e. exhibit exact 
area law.Q 

This reflects the striking behavior revealed by the numerical simulations 0, that 
over sufficiently large scales there is 'condensation' of vorticity carrying flux. Over a 
hypercube of side length of about 1 fm, the weighted probability at large j3 for nonzero 
(mod N) flux configurations goes to unity. The weighted probability then that one finds 
a vortex of at least this transverse size going through a 2-dim face on the hypercube 
boundary is approaching one. More generally, above this scale, vortices of any length, 
by corresponding appropriate adjustment in thickness, can occur at practically zero free 
energy cost. One may view this as percolation of vortices in the following sense. If one 
considers any two disjoint segments on the boundary of a large 4-dimensional simply 
connected region, the probability of being joined by a vortex of sufficient thickness (> 
1 fm) is finite. 

This picture of 'percolated' vortices in all possible [fiv] orientations, with flux in 
intersections being conserved only mod 2, implies that in general it is difficult to un- 
ambiguously identify individual vortices.0 Rather, in the absence of obstructions or 
boundaries introduced by external probes, one can talk about an average nonvanish- 
ing vorticity field defined on the coarse scale of 1 fm, measured by the 'circulation' 
(plaquettes, Wilson loops) above this scale. 

What simple effective theory can describe this vacuum? We stress that we mean an 
effective theory strictly of only these long range vortex fluctuations (confiners) resulting 
from integrating out to an appropriate scale. Even at that scale there will of course 
be all kinds of other SU(2) (SU(N)) fluctuations which we consider irrelevant for 
confinement. Let us list the minimal requirements for the effective theory: 

(i) On a coarse scale of about I = 1 — 1.2 fm, the partition function should be 
expressible solely in terms of vortex excitations (coclosed (closed dual) surfaces 
of codim 2). 

(ii) The mod 2 (mod N) property should be manifestly incorporated. 

7 Recall that it is a rigorous result that the potential cannot rise faster than linearly. 

8 The way embedded non-self-interacting 2-dimensional surfaces (in our case, surfaces of a certain 
thickness) can grow to densely fill 4-dimensional space, the so called capped gropes and towers, has 
been extensively studied in manifold theory pX[| . 
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(iii) The vortex flux through each coarse scale plaquette should incur an action equal 
to the vortex free energy per unit length for thickness ~ I, as defined and com- 
puted from (|3p. This amount of free energy for the confiners (their action on the 
coarser scale) is the cost above the sea of all other vacuum fluctuations. (Again, 
note that this depends only mod 2 on the number of 'vortex-introducing' singular 
gauge transformations injected in the box in (|3|).) 

(iv) On this coarse lattice, the thick vortices should be in a 'percolation phase'. Thus, 
despite (i), there should not be a useful expansion of the theory - i.e. a conver- 
gent, or, at least formally, systematic expansion scheme allowing computation of 
observables - such that each term in the expansion is characterized by a well- 
defined number of vortex excitations. 

(v) The Wilson loop expectation should give strict area law as in the bound (|TT|). 

Let then A c be the coarse lattice of lattice spacing ~ I, and Xp £ -^(2) variables 
residing on its plaquettes. Then the simplest way to implement (i) above is by the 
partition sum: 




The measure enforces the constraint 

II Xp = 1 ( 29 ) 

pddc 

on the plaquettes forming the boundary of every 3-dimensional cube c on A c , so only 
excitations on coclosed sets are allowed. Equivalently, on the dual lattice, (^) assumes 
the form 

n x; = i (so) 

p£d*b 

on the plaquettes forming the coboundary of every bond. The requirement (ii) is then 
automatically taken care of. 

The general form of the effective action A e ^ 

A eS = Pes J2x P + fc P Yl XpXp' + fop XpXp'Xp"^ (31) 

P (p,p')&d*b (p,p>,p")£d*b 

involves, in addition to the basic plaquette term, quasilocal interaction terms involving 
two or more plaquettes in the coboundary of each bond, etc. Now, from (iii), and (H), 
we must have 

/5 efr ~exp(-p(/5)/ 2 ) (32) 

giving, in principle, the coupling /3 e g in terms of the coupling at the original lattice 
spacing, as p must approach the string tension for sufficiently large /. From the numer- 
ical simulations |3, I ~ 1.1 fm. This gives (3 e g ~ 0.002. This very small value reflects 
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of course the fact that at this choice of the physical length I vortex flux is found to 
become very 'light'. Correspondingly, the terms involving products of two or more 
plaquettes must be of order /3^ s and higher, hence entirely negligible. 

The effective model fl28| ) is now seen to simply be a Z{2) LGT. Indeed, the constraint 



(E9|) can be explicitly solved by introducing Z(2) bond variables 75 by: 



X P = II 7b • (33) 

bedp 

Then 

Z ^= II d "1b exp (flea J2 Idp H ) , (34) 

JZ(2) Y P 

where the ellipses indicate the additional clover and higher loop terms corresponding to 



the additional terms in (pl|) . The theory is in the deep strong coupling regime /3 e ff "C 1. 
Thus (iv) above is indeed satisfied. The theory can be treated in the strong coupling 
expansion. It cannot, however, be meaningfully expanded in its vortex excitations - 
that would be appropriate in the weak coupling (3 e R 3> 1 regime in the form of the 
usual weak coupling expansion for discrete groups (dilute vortex gas) [[llj]. 

It should perhaps be explicitly pointed out that the Z{2) variables in fl2"8|), (p4[), 
serve as an effective description of long distance fluctuations creating topological Z{2) 
flux (elements of tti(SU(2)/Z(2)) in the original theory [|T^]. They have nothing to do 
with the Z{2) part of the original SU(2) bond variables. Note that the Z(2) gauge 
theory interaction, together with (^), correctly reproduce the effects of flux spreading 
and thickening of vortices while conserving flux mod 2. Correspondingly, the Wilson 
loop now automatically gives the correct area law. 

The Wilson loop in fl3~4|) represents the coupling of an external quark current to the 
long distance confining fluctuations. Its replacement by dynamical quarks introduces 
fermions in the medium of these fluctuations. Since the effective coupling is strong, it 
induces dynamical chiral symmetry breaking (CSB). 

CSB in strongly coupled LGT has actually been demonstrated analytically in the 
superstrong gauge coupling limit (no plaquette action term) by expansion or mean 
field (large N or d) approximations []13~| |, [|L4}| , and rigorously by infrared bounds |15| 



It is physically obvious that the result extends to a finite region in the strong coupling 
regime. [] 

We may try to use the results in fl3| - [|l^] to get an estimate of the contribution 
to the quark condensate in the effective theory. Corrections from the plaquette term 
in the action are totally negligible due to the smallness of /3 e fj. One then has: 

qq) = z(l)N i jf (l - k(d)) 1/2 , (35) 



9 It should be possible to prove this by cluster expansion techniques around the /3 — point, though 
not so straightforward for technical reasons (bounding terms with Grassmann integrands). 
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where § < fc(4) < 0.35. z(l) is some renormalization factor that, in a more sophisti- 
cated treatment, should depend on how fermions are introduced at the original lattice 
spacing. Here we naively set it equal to one - this is equivalent to simply taking stag- 
gered fermions on the coarse lattice. With I = 1.1 fm, this gives (qq^ = (195 MeV) 3 
for N = 2, and (223 MeV) 3 for N = 3. This indicates that the quark condensate may 
be entirely accounted for by the long range confining fluctuations. 

Following our previous development, the effective Z(N) theory appears to emerge 
rather naturally, and in fact in a fairly unique manner. The idea that an effective 
theory of long range vortex fluctuations must be a Z(N) LGT is not new, but has not, 
we believe, been formulated in this way before. Recently, a model equivalent to (^8|), 
in the representation ( PU[ ) and employing additive Z(2) variables, was considered in 
fl6fl , apparently without any reference to Z{2) LGT. 



5 Conclusions 



In the present paper we studied the energetics of how vortices can disorder Wilson 
loops of different sizes. Vortices of any thickness smaller than the linear size of a given 
Wilson loop can link with it and contribute to disordering its average. Here we pointed 
out that for Wilson loops of any given size it is the vortices of "maximal" thickness, 
i.e. thickness of the order of the linear size of the loop, that give the most important 
contribution resulting in an area-law suppression of large Wilson loops and a linear 
heavy quark potential. On the other hand, vortices of any fixed thickness contribute 
only with a logarithmically suppressed term to the potential. This is in contradiction 
with the naive picture of confinement by vortices which assumes that vortices of fixed 
thickness can link with a fixed probability with arbitrarily large Wilson loops. The 
correct picture must take properly into account the relation between the length and 
thickness of a vortex imposed by the free energy requirements for creating the vortex, 
as well as the interaction between vortices introduced by the constraint of mod N 
conservation of the vortex flux. This picture of vortices naturally yields a long distance 
effective Z(N) gauge theory above the confinement scale of around lfm. The effective 
theory is deep in the strong coupling regime which makes it impossible to interpret it in 
terms of a simple vortex gas expansion. The only useful expansion one can consider is 
the strong coupling one. Being deep in the strong coupling regime, the effective Z(N) 
gauge theory naturally produces chiral symmetry breaking in the presence of fermion 
fields. 
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